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ABSTRACT 

The spatial and temporal noise properties of diffuse radiation is investigated 
in the context of the cosmic microwave background (CMB), although generic 
formulae that enable application to any other forms of incoherent light of a 
prescribed energy spectrum are also provided. It is shown that the variance 
of fluctuations in the density and flux consists of two parts. First is a term from 
the spontaneous emission coefficient which is the contribution from a random 
gas of classical particles representing the corpuscular (photon) nature of light. 
Second is a term from the stimulated emission coefficient which leads to a 'wave 
noise', or more precisely the noise arising from the superposition of many plane 
waves of arbitrary phase - the normal modes of the radiation. The origin of this 
second term has never been elucidated before. We discussed one application. 
In the spatially homogeneous post-inflationary epochs when the universe was 
reheated to GUT temperatures, the thermal CMB density fluctuations on the 
horizon scale is of order the WMAP measured value of 10 -5 . Beyond (larger 
than) the horizon, the power spectrum of perturbations could moreover assume 
the observed form of P(k) ~ k if thermal diffusion to equilibrium does not merely 
involve classical particles, but also non-localized wave components. The outcome 
of our study clearly demonstrates that such components do indeed exist. 

1. Introduction: radiation noise and Poisson statistics 

The fluctuation properties of diffuse radiation is a topic of increasing importance, es- 
pecially with the advent of cosmological probes like COBE and WMAP, as they measure the 
density perturbations of the early universe that seeded the formation of structures (Bennett 
et al 2003, Hinshaw et al 2007, Spergel et al 2007). Yet, aside from the extrinsic physical 
processes that generated and grew these seeds, the cosmic microwave background (CMB), 
along with other forms of diffuse sky background, has its own intrinsic noise. It is helpful to 
acquire a better understanding of this noise, so that at the very least one has a firm handle 
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of what to expect even in the limit of a perfect instrument and a smooth sky (the latter from 
the 'extrinsic' standpoint). 

The motivation of our present paper is to take a close look at the intrinsic properties 
of radiation noise, with the black body radiation (BBR) as our prime focus, although gen- 
eralization will be made to include other types of diffuse emissions. In the field of optics, 
this investigation was initiated and vigorously pursued by several authors during the decade 
of 1950 - 60 (Purcell 1956, Hanbury Brown and Twiss 1957, Mandel 1958), and in recent 
times the work was continued mostly by researchers of ultrafast optics (Trebino 1997). Our 
purpose here is to further develop and articulate the findings of these authors in the context 
of astrophysical observations and interpretations. 

We begin by revisiting BBR itself, by considering lengthscales and timescales short 
enough that any effects due to the Hubble expansion can be neglected. The probability of 
having the integer n as occupation number of some normal mode j is governed by Bose- 
Einstein (BE) statistics, as 

P n = (1 - e-^) e~ n ^ , (1) 
where f3 = 1/(/cbT). The mean occupation number and its variance are, as a result, 

™i = e ffl^_i ' and ( 6n ^ = e ^-l ^ 



Similar parameters for the energy contribution from this mode are €j = rijhujj and (5e 

i 3 fh 2 u 2 



(5rij) 2 h 2 uj 2 . Since modes within any black body cavity of volume V are independent, the 



variance of the total energy e in this volume is simply (5e) 2 = J2j(,$ e j) 2 - 

Now, taking the continuum limit and applying the usual phase space density formula 
for the number of electromagnetic wave modes per interval dui of angular frequency 

r-2x 7 ^-/d 3 k, or%^=2x^-, (3) 

3 

where the factor of two represents the polarization degrees of freedom, one obtains 

87TV{k B T) 4 r° X 3 dx _87T 5 V{k B T) 4 

6 ~ {hcf J 1^1 ~ 15(/ic) 3 ' 

and 



(5ef 



8nV(k B T) 5 f°° x 4 e x dx 32ir 5 k B T 5 V 

J {e* - l) 2 ~ 15/i 3 c 3 ' 5 ' ' 



(he) 3 

Hence the contrast in e, or the energy density u, is 
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in agreement with the derivation of Peebles 1993 (S = 4e/(3T) is the entropy of the radia- 
tion). 

An important point of physics has to be made about eq. (J2J). The variance {Srij) 2 in rij is 
greater than the mean Hj, i.e. BBR noise does not precisely follow Poisson statistics. Yet the 
common misconception that in the complete (quantum) theory radiation noise arises solely 
from photon counting may still be defended. Thus, when integrated over all frequencies, 
eqs. dl]) and ([SD give the distinctly Poissonian behavior of (Se/e) 2 ~ 1/N where N = J dn is 
the total number of photons in the cavity of volume V. Also, for the high frequency modes 
x = fkj/(kBT) ^> 1, (5n) 2 pa n, i.e. the discrete nature of radiation as randomly moving 
photons prevails there. However, by writing the variance as two terms3, 

(5n) 2 = n(n + 1), or (5n) 2 = j~^ 2 = + J^T^r ( 6 ) 

we see that the Poisson noise term n dominates at x ^> 1, whereas the n 2 (rightmost) 
term is significant only at low frequencies x -C 1. Focussing upon the latter, we may 
apply the x <Cl limit to eq. (j3J) and obtain, for a frequency interval 5x in this regime, 
(5e) 2 ~ x 2 5x/h 3 ~ ^Sv/^bT) 3 , so that the variance is no longer dependent upon h. Clearly 
here, we are dealing with noise of a 'classical' nature. But what exactly is it? The question 
has never been answered rigorously. 

2. The precise origin of the non-Poisson component and its importance 

One could take the questions raised above as a red herring, or at best pedagogical: what 
prevents one from simply appealing to radiation as an ensemble of photons, which in turn 
are bosons with the manifestly non-Poisson distribution of eq. flTJ? In other words, one 
already has an exact theory of BE statistics that renders any 'analysis' like the previous 
paragraph quite obsolete. The answer is three-fold. Firstly, the application of BE statistics 
to systems of various sizes is non-trivial even for BBR. Results like eqs. (j3J) and (jSJ) actually 
describe an incoherent ensemble of ~ V/(AX) 3 'pure BE' systems with AA denoting the 
bandwidth of the spectrum (i.e. the range of wavelengths within which most of the emitted 
flux resides); each of such systems corresponds to an elementary cell in the standard phase 
space quantization scheme. If e.g. the volume of space in question is < A 3 , coherence effects 
will play a role, modifying the variance from the form eqs. (H| or (131), see Mandel 1958 
(also Purcell 1956; Hanbury Brown, and Twiss 1957). The same 'partitioning' phenomenon 
applies as well to intensity fluctuations 61 over time intervals large or small compared with 
1/Au. Such behavior deserves further investigation, which is one purpose of this paper. 
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Secondly, from the standpoint of modeling cosmological structures there is especially 
a need to press for better understanding of the two terms in eq. ([6]). The primordial 
density perturbations that seeded structure formation have hitherto been characterized by 
the two-point correlation function £(i"i — r 2 ), the Fourier transform of which is the power 
spectrum P(k), which describes how particles are positioned in space via a locally random 
(Poisson) process. If, however, during the radiation era the fluctuations have a component 
that cannot be described in terms of the distribution of classical Poisson particles at all, 
then the properties of £(ri — r 2 ) must be inferred without this perspective, i.e. an altogether 
different microscopic treatment will be necessary instead. 

Thirdly, in astrophysics and experimental optics one could easily envisage other types of 
radiation, such as the diffuse X-ray background, or simply any arriving signals of incoherent 
light - often known as white light. If the mean energy spectrum (more precisely spectral 
energy density) U(u) of the radiation is known, can we infer the fluctuations in the energy 
density at any length scale, or the intensity measured over any time scale? Obviously, a clear 
picture of the origin of all the noises is prerequisite to the solution. 

To elaborate upon this third point, let the mean occupation number per mode be n{oS), 
so the mean energy density of all modes is 

1 



U 



7T 2 C 3 



ui hun(uj)du! 



U(uj)doj, 



(7) 



where, from eq. (jSJ) 



UJ" 



7T 2 C 3 
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Now the distribution n(uj) assumes, for a given value n(u), the most likely occurring form 
when its Entropy is maximized subject to this constraint on the mean, irrespective of whether 
each mode is in 'internal equilibrium' or not (if it is, the 'most likely' qualifier will be 
unnecessary, see Jaynes 1957 for the complete theory). Under the scenario of states being 
occupied with complete a priori freedom, i.e. boson systems with no prior restriction on n, 
the resulting n(u) would have the variance of [5n(cj)] 2 = n(u)[n(u) + 1]. Hence, as in eq. 
©, 

1 f°° 1 
(SU) 2 = ^ r — LU 2 (Huj) 2 n(u)[n(u) + l)}duj 



n 2 c 3y 



If the radiation is BBR, i.e. 
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hwU(uj) 
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U{u) = U{uj) 



ftjjj" 



1 



duo. (9) 



(10) 



tt 2 c 3 e^-V 

it will be trivial to show that eq. (jUJ) is then the same as eq. ( 
non-Poisson (or 'wave') noise given by the U(uj) integral and the U 2 {uj) integral respectively. 



with the Poisson and 



- 5 - 



In spite of the elegance of the above derivation, one still ought to seek a more substantive 
physical justification of eq. Qj that does not appeal to any line of reasoning based upon the 
'game of chances' of Jaynes 1957. We only resorted to such an argument because of our lack 
of fundamental understanding of the non-Poisson noise. 

3. A heuristic treatment of wave superposition noise 

The 'n 2 -noise' contribution to the total variance of BBR fluctuations may further be 
elucidated, even if still heuristically, by Maxwell's wave model of the radiation. Let BBR 
modes exist in the usual manner within some volume V. The normalized electric field of 
mode j is written as 



where ej is the polarization unit vector. The total energy is e = uV = VY2j \Ej\ 2 - It seems 
reasonable to assume that Ej are uncorrelated random variables, with a gaussian distribution 
of zero mean and variance (\Ej\ 2 ) ~ (ej)/V, where (ej) is as defined in the text following eq. 
(J2J). Thus, the variance (5e) 2 is the sum of the individual (Sej) 2 , i.e. 



3 

For a complex gaussian random variable of zero mean, the second last term is twice the last 
term, so that 



which is the integrated contribution to the total variance from the 'non Poisson' (rightmost) 
term of eq. (jSJ). 

It should be clarified that classically the random complex nature of the field amplitude 
Ej is due to the wave phase <j>j (which changes from mode to mode without correlation) 
rather than the photon counting noise, which is why the variance is not associated with 
Poisson statistics. In fact, we can rewrite eq. ([7j) as 



Ej(r, t) = EjBj exp i(kj ■ r — uijt); E 3 -(r, t = 0) = Ejej exp(zk,- • r) 



(11) 




(12) 




(13) 




(14) 



where \Ej\ 2 = U(uj) and U is as defined in eq. f fTUj) . In the next section we proceed to 
perform a more rigorous calculation of this wave noise (or, better still, phase noise), by first 
attempting just the one dimensional problem of intensity fluctuations at various time scales. 
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4. Intensity fluctuation due to random mode phases, the coherence length of 
BBR 

The total intensity (or flux, in ergs cm -2 s^ 1 ) of, for simplicity, a collimated beam of 
BBR arriving at the observer's point r = is given by 

i r /2 

I(t) = - / c\E{t)\ 2 dt, (15) 

where r is the time interval over which the intensity is measured, and 

E(*) = ^E i (t), (16) 

3 

withB Ej(t) = Ej(r = 0,t) and the right side is as defined in eq. (fT4j) . 

We shall soon see that J(r) is a sum of two terms. One is a constant 'DC term, 
independent of r in particular, while the other is a random 'AC term of zero mean, the 
variance of which decreases with increasing r. This second term causes J(r) as evaluated 
by means of eq. ( JTBT) with different 'time origins' to fluctuate from one interval of duration 
t to the next. The behavior is actually well known in optics (Trebino et al 1997), and is a 
consequence of the Uncertainty Principle. To see this readily, one could return to eq. (|T4l) 
and set 

<f>j = for all j. (17) 

The sum over j performed in eq. (TT6T) to obtain the total wave electric field would then 
be, in essence, a Fourier inverse transform of the Planck function (more precisely v/Z7(o>)). 
The result is a 'pulse' of duration At ~ 1/Au, where Auj is the spectral width of U(u), as 
is apparent when one plots the intensity profile I(r) against r with high time resolution, 
Figure 1. This scenario is known in the field of Fourier optics as a 'transform limited pulse', 
viz. the Uncertainty Principle is satisfied in the 'minimal wave packet' sense of Auj At ~ 1, 
and it is attained by setting all the mode phases to the same value - the highest degree of 
spectral phase coherence that can only be found in laser light. 

If, as is normally the case of ordinary 'white light', the normal modes have more com- 
plicated phase relationships than eq. (fl~7]) . the corresponding pulse will in general last longer 
than the minimum duration, i.e. it will satisfy AuAt > 1. When the other end of the 
extreme is reached, such that the phases <pj exhibit the maximum complexity by assuming 



2 We shall see in section 6 that the validity of choosing a hxed point r = to do the calculation requires the 
detector to be centered at this point, and of size < one radiation wavelength (to maintain spatial coherence) 
which in the case of 3K light is ~ 3 cm. 
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uncorrelated random numbers as their values, one is dealing with 'ideal white light' (BBR 
falls under this category), viz. an arbitrarily long and continuous puls^f] with AcuAt ^> 1 
(Trebino et al 1997). Nevertheless, this continuous light train is in detail broken into many 
'spikes' of all widths, ranging from one period of electric field oscillation at the BBR peak 
frequency to the full span of L/c. Such 'spikes' are to be realized in the 'AC term of the 
previous paragraph. 

To substantiate these points, we numerically computed eqs. (1T11) through (TTB]) at two 
different resolutions r, one beneath the reciprocal bandwidth of the 3K BBR spectrum (i.e. 
t <C 1/Au) and the other well above it, this time assigning random phases <pj to each mode 
j as is appropriate for the BBR. The light intensity averaged over temporally contiguous 
intervals^ 4 ] of r, I(r), is indeed now seen in Figure lc and Figure 2 to be continuous, yet it is 
equally obvious that there are hardly any changes in / on the scale of t < 10~ 13 s; when r 
is larger, at e.g. r = 1CT 10 s, noises then become apparent. In fact, as will be proved below, 
variations occur on all scales inside the range 

7^<r<^ (18) 
Au olo 

and with properties determined solely by the BBR spectrum B(u) and independent of the 
mode quantization frequency 5u used in our simulation. For values of r approaching the left 
inequality, the light curve becomes relaxed and regular, eventually altogether steady when 
t < 1/Auj as Figure 2a reveals. This is consistent (again) with the Uncertainty Principle, 
viz. on timescales At < 1/Au there can be no randomness in the signal. Far from this 
limit, the intensity turns erratic as the noise decorrelates from one resolution element r to 
the nexl0, Figure 2b. It should also be mentioned that the right inequality of eq. ({TBI) , viz. 
8uj < 1/r, is equivalent to the constraint imposed by the periodic boundary condition as 
applied to a cube of side L = ct. 

Our simulations for Figures 1 and 2 were pursued with a mode quantization Su = 10 10 
rad/s. Further beyond plotting the intensity at the two r values discussed, we calculated 
from such plots the normalized variance (SI /I) 2 (where I is the 'global' mean intensity which 
is a constant independent of r) at each scale between r = 10" 14 s and r = 10~ 10 s. The 



3 This is illustrated in Figure 11 of http://en.wikipedia.org/wiki/Coherence_(physics) . 

4 To obtain I(t) for many of these intervals of r we need to repeatedly compute eq. (|15[) using a fixed set 
of random phases but shifting the time of the mid-point of the integration limits by the amount r for every 
interval, alternatively one can keep re-randomizing the phases on a per interval basis. 

5 Figure 2b bears resemblance to the fluctuation pattern of 'white light' from incoherent mode superposi- 
tion as depicted in http://en.wikipedia.Org/wiki/Image:Spectral_coherence_continuous.png. 
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Fig. 1. — Two simulated intensity profiles (arbitrary units and photon Poisson noise ex- 
cluded) of unpolarized light with the 3K black body spectrum of la (extreme left). The 
detector is assumed to be < 1 cm, or one radiation wavelength, in width and length, to 
avoid complications due to spatial incoherence. Constituent normal modes are spaced 5u = 
10 10 rad/s apart, lb (middle): all the modes are in phase, resulting in a single pulse which 
is the Fourier transform of the spectrum. The spectral width Auj of la and the pulse width 
At of lb satisfy Auj At ~ 1. This is the scenario of a 'transform limited pulse' (or 'minimal 
wave packet'), and occurs only when the modes have full phase correlation, lc (right): the 
modes are randomly phased without any correlation among them; here we see an enduring 
light curve with aperiodic spikes and structures on all scales. Both profiles were constructed 
using eqs. (fT^|) through (TL6]) . with a timing resolution of r = 10~ 14 s. 

results are shown in Figure 3 as the data points of a power spectrum (in the WMAP sense of 
the words, i.e. fluctuation power at various scales). Towards r = 1CT 10 s the variance ~ 1/r, 
suggesting (as already said) that neighboring bins of r decorrelate when r becomes large. 
The argument is as follows. From eq. ( [15]) the average intensity over some long interval 
T = Nt may be written as 

1 



it = c x — 



/ \E(t)\ 2 dt+ / \E(t)\ 2 dt + ■ ■ ■ + / \B(t)\ 2 dt . (19) 

Jo Jt J(N-1)t 

If each of the N terms on the right side fluctuates freely with a variance a 2 ., we will have 
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Fig. 2. — Intensity fluctuations (phase noise component) of 3K BBR as seen in two resolu- 
tions, 1CT 13 s and 1CT 10 s. Mode frequency spacing adopted by our simulation is 6u = 10 10 
rad/s as in Figure 1. 

(SI r y = Na 2 /T 2 = cr 2 /(Nr 2 ); but since (5I T ) 2 = a 2 /r 2 we obtain (5I T ) 2 = (5I T ) 2 /N, or 
T(5Ir) 2 = t(5I t ) 2 = constant. Hence (Sir) 2 ~ This trend holds as long as we compare 
two timescales in which the noise in the larger bin is simply the incoherent reinforcement of 
those in the smaller ones. 

At r < 10~ n s, however, we see in Figure 3 a transition towards different behavior. 
As Figures lc and 2a showed, the intensity changes over short intervals are regular and 
highly correlated, and variances no longer add. Eventually, (Sir) 2 = (SI T ) 2 and the power 
spectrum flattens at the maximum, i.e. no 'extra' fluctuations may be found by measuring 
the intensity using more and more accurate clocks, so 51 does not rise any further. From 
Figure 3, one finds that the cross over point from constant 5I T to the 8I T ~ 1/\/t dependence 
is at r = t c pa 1CT 10 s. This length may be used to define the coherence scale of BBR at 
T = 3 K. Beneath r c (or cr c = 3 cm) a wavetrain of BBR exhibits good phase correlation, 
above it the phase becomes randomized. In another manner of expression, 3K radiation is 
capable of maintaining the integrity of its electromagnetic oscillation pattern for u; max r c ps 
100 cycles. The overall picture is that one could divide space-time into many cells of this 
critical dimension, and assume coherence within each cell. Thus e.g. the point raised earlier 
in footnote 2 may now be understood in better terms: the power spectrum simulated and 
derived in this section is valid only for a detector size of < 3 cm because we have not taken 
into account the additional effects of spatial incoherence. 
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time scale (s) 



Fig. 3. — Power spectrum of temporal intensity fluctuations of 3K BBR, excluding photon 
Poisson noise. The uncertainties in the numerical simulations are < the size of the circles, 
and are determined by generating data from 350 different sets of random phases (f)j, then 
evaluating the scatter among the variances that result from each. The mode quantization 
adopted throughout is 5uj = 10 10 rad/s. 

Finally, it is also clear that the salient features of this section applies to diffuse radiation 
of any spectrum B(u). 

5. Analytical proof of temporal phase noise variance 

The outcome of section 4 is not only relevant to BBR, but to incoherent light (or 'white 
light') with any broadband spectrum. A more rigorous mathematical proof to support the 
simulations is clearly desirable. 

The intensity J(r) at resolution r, eq. (fT5]) . has a 'DC and an 'AC contribution as 
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mentioned after eq. (fT6j) . These may now be written as 

/ DC = Vfi J ; J A c(t) = V Y] e kitn -e k /, n V5 fe B fc /-4 i-J-e^.m,,, (20) 

i „ , TT^i \ (uJ k -uJk')T 

] k,m (k',n)^(k,m) v ' 

where </>k,k',m,n = 0k,m — 0k',n ( m , n = 1,2 are polarization labels), B k = cU k = hkc 2 / {ef 3 *** — 
1). Evidently Jdc is just a constant. Pairing each term in the series of Iac with its complex 
conjugate, the factor in eq. (f2"0|) is then replaced by 2cos0, and the J^ k , procedure is to 
be done with the understanding that only one member of every (k, k') and (k', k) pair should 
appear in the sum. The mean of Jac is evidently zero. The variance is, in the continuum 
representation where B(u) = cU(u) with the latter given by eq. (jlOp . 



2 

9 r°° f°° sm - — 
(6I T ) 2 = - / dujB(uj) / duj'B(uj')—± rjJ-, (21) 

T Jo Jo — 

after applying a correction factor of 1/2 to free the step from the awkward 'permutation' 
restriction above; also the relation (cos 2 </>) = 1/2 was used for the random phase differences 
0, and the sum over random polarizations follows the standard formula 

E^(e k , m -e k , n ) 2 = l+(^) =2 (22) 

(k-k ; « kk' because u' ~ c<j as we saw, and k||k' for a collimated beam), this sum replaces the 
factor of four from the helicity enumeration in the two phase space factors that accompany 
B k and B k > as they become B(u) and B(u'). 

A plot of (6I T /I) 2 against r, using eq. ff2T]) and the 3K Planck function of eq. ffTUj) . is 
given by the solid line of Figure 3, where excellent agreement with data from the numerical 
simulation of section 3 is apparent. There are two obvious limiting cases in which eq. (12"T1) 
simplifies. First is when r 1/Au, the reciprocal of the BBR spectral width. In the second 
integral we may now write sin(a/ — uS)r ~ (u/ — uj)t, and {SI T ) 2 becomes a constant, as can 
be seen from the flattened top portion of the variance curve in Figure 3. Second is when 
r 3> l/Au. Here it is tempting to use the property sin 2 -?/' < 1 to recast eq. fl2"T|) as 

(5I T ) 2 < - / duB(u) / dcu' / >— (23) 



(uj 1 — uj)t 



r 2 Jo ' Jo ' (w' -uj) 2 



to conclude that (SI T ) 2 must asymptotically fall with r at least as steeply as 1/t 2 . This is 
wrong for the following reason. The function B(u') / (u/ — uj) 2 always diverges at some point 
in the uj' integration in such a way as to yield an infinite upper limit to (5I T ) 2 , of the form 
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Instead, to arrive at the correct asymptote one must observe that for r ^> 1/Au the 
J du' step in eq. ( I2ip converges rapidly: within some interval u = \u' — uj\ <C Au, enabling 
us to write B(oj') m B(u) and take the function B(u) out of the u' integration, i.e. eq. (12"T1) 
reads 



w ° sin 2 (^f) 2 f°° 2 f°°sm 2 f 1 



(5I T ) 2 = 2 / du[B(u)} 2 / du ^ l 2 > = - / du[B(u)] 2 / — ^df oc -, (24) 
JO J -wo u T T Jo Jo 4 r 

for t ^> 1/ Au, where u = u' — u and £ = ujt/2. Note that the integral over £ was taken 
to £ = oo, resulting in the pure dimensionless number of tt/2, because when r is large 
convergence to this number only requires an excursion of u across a very small passband 
(~ Uq in width) within the Planck spectrum. The asymptotic 1/r dependence of the variance 
is, finally, 

f5I T \ 2 225h f°° x 6 dx 9 h / r x^ 1 . 

( —1 = — / = 7.3 x 10 2 x - — = 0.0123 ( -— I for T = 3 K, 

\I J 27r 8 rkT J (e^-l) 2 7r 8 rkT VlO- 10 sJ 

in agreement with the simulated data of Figure 3, and consistent with the rationale presented 
around eq. (1191) of large scale fluctuations being fed by incoherent noises at smaller scales 
(see also Mandel 1958). 

It is worth a mention in passing that one can also pursue the simpler problem of the 
variance in E(t), where 

1 r r/2 ^ sin 
E(t) = - / Re[E(t)]dt = V2v/5~ ^^cos^- (25) 

is the magnitude of the BBR electric field averaged over the time r (a quantity that obviously 
fluctuates about zero mean), viz. 

mr = j22 B ;-^ffi, (26) 

U A T 
3 3 

with the last result obtained by applying (cos 2 <p) = 1/2 as before. The main difference from 
the {SI r ) 2 calculation of eqs. ( 12~TT) through (1231 is that this time, we can actually follow the 
approach of eq. (1231) without making a mistake, viz. by writing eq. fl26|) as 

mr<^ 2 r^^M. (27) 

r 2 Jo u 2 



The key point is that there is no longer any infra-red divergence in the integrand: at low 
frequencies B(u)/u 2 remains finite due to the Rayleigh- Jeans Law B(u) ~ u 2 for small u. 
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The large scale r.m.s. electric field, then, can be rationalized in terms of the arguments 
presented near eq. (123]) . i.e. it falls off at the rate of 1/r 2 or steeper, and the velocity of any 
charged particle v = (eE rms /m)r does not increase with r. This is a reflection of the thermal 
equilibrium environment of BBR: the long wavelength fields have already been dissipated by 
re-absorption and ensuing particle acceleration. 



6. Spatial contrast of the energy density of BBR - the phase noise contribution 

We are now in a position to tackle the problem of spatial fluctuations - the task identified 
in sections 1 and 2 of finding the precise origin of the phase noise contribution to the spatial 
variation of BBR energy density. This time, we set in eq. (1141) t = 0, i.e. we take a snapshot of 
the radiation filled space (more precisely defined as an interval of duration < 1/Au centered 
at the time origin) and examine the density distribution everywhere. The energy density 
U(V) measured over some cubical volume of side L is given by 



«(V) = ~ J |E(r)| 2 rf 3 r, {2i 



where V = L 3 . Following the notation of sections 4 and 5, this means u{V) = udc(V) + 
u\c(y), where 

ubc(V) = U(k); u AC (V) = ^J2 E e kim -e kl , n ^U(k)U { k') j e^~ k '^d 3 r, 

k,m k.m (k,m)^(k',n)^(6fc,m) 

(29) 

with <p = 4>(k, m) — <p(k', n), m,n = 1,2 are polarization labels, and 

& w = < 30 > 

(U is the same as U of eq. (TTOj) with the phase space factor removed to enable a full three 
dimensional mode sum). 

The integral in eq. (T2"9"|) may then be performed. After pairing with conjugate terms, 
we obtain 

«ac(^) = 16> > e k , m • ev in JU(k)U(k') f — cos0, (31) 

tig-Lj K/ylj Kj Z ±J 

k,m (k',n)^(k,m) ' 

where k, = k' — k, and it is understood that for every pair of terms differing from each 
other by a mere swap of mode ordering, one member of the pair drops out of the sum. The 
variance is 

{5uf = 256^ U{k)U(k') ((e k , m • e k / in ) 2 )(cos 2 < 

k,m (k',n)^(k,m) 



sm(n x L) sm{K y L) sm(n z L) 

Kj^L KjyL hc^L 



(32) 
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Next, in eq. (1321) the sum over k' is replaced by J <i 3 k' = J d 3 K, = J dK x dK y dK z , and the 
resulting integration of the three sync functions can be done by means of the reasoning 
that led to eq. (1241) . i.e. for large L only the range of k near n = (or k' « k) contributes. 
The outcome is, 



(8u) 



7T 2 C 3 



V 



POO 

/ do; 




Jo 





(33) 



where U(w) is given by eq. (flQj) and the polarization sum in eq. (132]) is again carried out 
means of eq. (122|) . i.e. since k « k' the double summation of eq. (122]) yields two as the 
answer, so that the 'helicity factor' of four in U 2 {uj) is reduced to two. The integral of eq. 
(13"3"|) may be recast as 

, c , 2 8vrA; 5 T 5 Z" 00 x 4 dx 
(5w) 2 



h 3 c 3 V J (e* - l) 2 

Referring back to eqs. (j3j) and ([H]), and bearing in mind that e = it/V, we proved that the 
'non Poisson' (rightmost) term of eq. ([6]), when integrated over all frequencies x according 
to eq. PI, agrees exactly with the phase noise variance of the energy density as calculated 
here. 

Further, it should be emphasized that the function U (u) in eq. (I3"3"j) can be the spectral 
energy density lA{u) of any diffuse radiation, and the equation will still yield the correct 
phase noise variance. Moreover, it is also evident that the Poisson noise (second last) term 
of eq. ([B]), has the variance of 

1 f°° 

(Su) 2 = - HuU(u)duj, (34) 
* Jo 

and the total variance is the sum of the two contributions, in agreement with the heuristic 
(maximum Entropy) treatment of eq. (Q. 

7. Conclusion: cosmological implications 

Turning finally to the CMB and cosmology, the observed CMB anisotropy over length- 
scales of Mpc or above is far in excess of the intrinsic 3K BBR noise discussed throughout 
this paper. It is widely believed that the origin of the CMB anisotropy is vacuum fluctuations 
generated during inflation, and caution must then be taken when one contemplates those 
very early epochs. Thus, as explained in Lieu & Kibble (2009), the intrinsic BBR noise at 
the end of inflation when the universe was reheated to kT ~ 10 15 eV (GUT) temperatures 
is by no means negligible; in fact, on the scale of the GUT horizon the density contrast 8e/e 
due solely to this thermal noise (assuming each island universe is causally connected within 
its own Hubble horizon and has therefore the time to reach thermodynamic equilibrium) is 
~ 10~ 5 , on par with the observed value of the same for horizon crossing modes of density 



perturbation at relatively much lower redshifts. This point is particularly significant since 
5e/e on the horizon is, for the observed Harrison-Zel'dovich form of the power spectrum 
P(k) ~ k of superhorizon fluctuations (for the meaning of P{k) see section 2), both a gauge 
and epoch independent quantity having a value set by the initial condition after reheating. 

Apart from the above justification, which concerns the amplitude of the density seeds, 
the furthering of our insights on the BBR noise as developed in this paper also helps us to 
better define the role of thermal diffusion. As case in point, one important application is 
to the problem of how P{k) evolves from an initially spatially homogeneous configuration 
(the highly ordered state that ensued from the reheating of a post-inflationary universe) 
as a result of the gradual thermalization of superhorizon regions. It has been known for 
some time (Zel'dovich 1965, Peebles 1974, Zel'dovich and Novikov 1983) that such a process 
gives rise to the power spectrum P(k) ~ k n where n > 4. The inequality for n is often 
referred to as the Zel'dovich bound. However, a careful examination (Gabrielli et al 2004 
and references therein) reveals that the underlying assumption adopted by these historical 
authors is the random walk of classical Poisson-like particles. Moreover, Gabrielli et al 
pointed out that if the diffusion mechanism at work in the very early universe involves 
also the 'fuzzy' wave-like properties of its constituent 'particles' - properties that 'blurs' the 
sharpness of the cutoff of interactions beyond the horizon separation distance because the 
'particles' cannot be localized as readily - the Zel'dovich bound might be broken, opening 
the range of allowed n to even include the Harrison-Zel'dovich value of n pa 1. Thus, in this 
context the present paper clearly plays a role in bolstering the argument of Gabrielli et al, 
viz. a radiation dominated ensemble does exhibit density fluctuations that originate from 
the phase interference of waves, hence it is by no means a system of classically diffusing 
particles. As explained in Lieu & Kibble (2009), under the circumstance the attractive 
prospect of accounting for both the amplitude and index of primordial density perturbations 
without invoking fluctuations in the vacuum energy then avails itself. 
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